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Abstract. I report recent results of numerical simulations designed 
to study the inner structure of Cold Dark Matter (CDM) halos. This 
work confirms the proposal of Navarro, Frenk & White (NFW) that the 
shape of ACDM halo mass profiles differs strongly from a power law 
and is approximately independent of mass. The logarithmic slope of the 
spherically-averaged density profile, as measured by (3 = — din p/dlnr, 
decreases monotonically towards the center, becomes shallower than isother- 
mal (/3 < 2) inside a characteristic radius, r_2, and shows no evidence for 
convergence to a well-defined asymptotic value (/?o) at the center. The 
dark mass contained within the innermost radius resolved by the simula- 
tions places strong constraints on /3q; cusps as steep as r -1 - 5 are clearly 
ruled out in our highest resolution simulations. A profile where the radial 
dependence of the slope is a simple power law, /3(r) oc r a , approximates 
the structure of halos better than the NFW profile and so may minimize 
errors when extrapolating our results inward to radii not yet reliably 
probed by numerical simulations. We compare the spherically-averaged 
circular velocity (V c ) profiles to the rotation curves of low surface bright- 
ness (LSB) galaxies in the samples of de Blok et al and Swaters et al. 
The V c profiles of simulated CDM halos are generally consistent with the 
rotation curves of LSB galaxies, but there are also some clearly discrepant 
cases. This disagreement has been interpreted as excluding the presence 
of cusps, but it may also just reflect the difference between circular ve- 
locity and rotation speed likely to arise in gaseous disks embedded within 
realistic triaxial halos. 



1. Introduction 

The pioneering efforts of Frenk et al (1988), Dubinski and Carlberg (1991), 
and Crone et al (1994), among others, led to the identification of a number of 
key features of the structure of dark matter halos assembled through hierarchical 
clustering. One important result of this early work concerns the remarkable sim- 
ilarity ("universality") in the structure of dark matter halos of widely different 
mass. This was first proposed by Navarro, Frenk & White (1996, 1997; here- 
after NFW), who suggested a simple fitting formula to describe the spherically- 
averaged density profiles of dark matter halos, 

P(r) = (r/r.)(l + (r/r a ))2- (1) 
1 
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A second result concerns the absence of a well defined central "core" of constant 
density in virialized CDM halos: the dark matter density grows apparently 
without bounds toward the center of the halo. 

Subsequent work has generally confirmed these trends, but has also high- 
lighted potentially important deviations from the NFW fitting formula. In par- 
ticular, Fukushige Sz Makino (2001), as well as Moore and collaborators (Moore 
et al 1998, 1999), have reported that NFW fits to their simulated halos (of much 
higher mass and spatial resolution than the original NFW work) underestimate 
the dark matter density in the inner regions (r < r s ). These authors interpreted 
the disagreement as indicative of inner density "cusps" steeper than the NFW 
profile and advocated a simple modification to the NFW formula (referred to 
hereafter as the M99 profile) with asymptotic central slope (3q = 1.5. 

The actual value of (3q is still being hotly debated in the literature (Jing 
& Suto 2000, Klypin et al 2001, Taylor & Navarro 2001, Power et al 2003, 
Fukushige, Kawai & Makino 2003, Hayashi et al 2003, Navarro et al 2003), 
but there is general consensus that CDM halos are indeed "cuspy". This has 
been recognized as an important result, since the rotation curves of many disk 
galaxies, and in particular of low surface brightness (LSB) systems, appear to 
indicate the presence of an extended region of constant dark matter density: a 
dark matter "core" (Flores & Primack 1994, Moore 1994). 

Unfortunately, rotation curve constraints are strongest just where numeri- 
cal simulations are least reliable. Resolving CDM halos down to the kpc scales 
probed by the innermost points of rotation curves poses a significant compu- 
tational challenge that has been met in very few of the simulations published 
to date. Indeed, rotation curves are generally compared with extrapolations of 
the simulation data that rely heavily on the applicability and accuracy of fit- 
ting formulae such as the NFW profile to regions that may be compromised by 
numerical artifact. 

The theoretical debate on the asymptotic central slope of the dark matter 
density profile, /3o, has also led at times to unwarranted emphasis on the very 
inner region of the rotation curve datasets, rather than on an proper appraisal of 
the data over its full radial extent. De Block et al (2001), for example, attempt to 
derive constraints on (3q from the innermost few points of their rotation curves, 
and conclude that /3q ~ for most galaxies in their sample. However, their 
analysis focuses on the regions most severely affected by non-circular motions, 
seeing, misalignments and slit offsets, and other effects that limit the accuracy of 
circular velocity estimates. It is perhaps not surprising, then, that other studies 
(van den Bosch et al 2000, Swaters et al 2003) have disputed the conclusiveness 
of these findings. The disagreement is compounded by the results of the latest 
cosmological N-body simulations (Power et al 2003, Hayashi et al 2003, Navarro 
et al 2003), which find scant evidence for a well defined value of (5q in simulated 
CDM halos. Given these difficulties, focusing the theoretical or observational 
analysis on (3q seems unwise. We shall try here to improve upon previous work 
by comparing directly the results of our simulations with the full radial extent 
of the rotation curves of LSB galaxies. 

I report here results from a major computational effort pursued by an in- 
ternational collaboration bringing together computational cosmology groups in 
Victoria (Canada), Seattle (USA), Durham (UK), and Garching (Germany). 
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Figure 1. (a) Density profiles of 19 ACDM halos, scaled to the ra- 
dius where the logarithmic slope takes the "isothermal" value of j3 = 
—dlogp/dlogr = 2. Thick solid curves show the NFW profile; dotted curves 
correspond to the M99 profile, (b) Circular velocity profiles, scaled to the peak 
of the V c curve. Note that the NFW and M99 profile appear to "bracket" the 
most extreme mass profile shapes of ACDM halos. 



We have now completed a series of about 20 simulations of ACDM halos span- 
ning roughly five decades in mass, from dwarf galaxies with circular velocities 
V c ~ 50 km s™ 1 to galaxy cluster systems with V c ~ 1500 km s _1 . This series 
follows a thorough convergence study where the role of all relevant numerical 
parameters has been assessed in detail (Power et al 2003). As a result, we are 
able to identify in each simulated halo a minimum convergence radius (r conv ) 
beyond which the mass profile is reliably reproduced: circular velocity estimates 
are accurate to better than 10% for r > r conv . Each of the simulated halos has 
several million particles within the virial radius, allowing for the mass profiles 
to be traced reliably down to within 1% of the virial radius. 



2. CDM Halo Density Profiles: Universality and Cusp Structure 

Figure 1 shows the density and circular velocity profiles of nineteen ACDM 
halos simulated in our series. The "universality" in the density profile shape is 
displayed by scaling each profile to r_2, the radius where the logarithmic slope 
is /3 = —dlogp/dlogr = 2. The left panel of Figure 1 shows that there is little 
difference in the scaled density profiles of CDM halos differing by up to five 
decades in mass. The same applies to the V c profiles, which in the right panel of 
Figure 1 have been scaled to the peak in the circular velocity curve. The fitting 
formulae proposed by NFW or M99 are fixed curves in these scaled units, and 
are seen to more or less "bracket" the most extreme shapes of the mass profiles 
of the simulated halos. 
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Figure 2. (a) Logarithmic slope of the density profile as a function of 
radius for (from left to right) "dwarf", "galaxy-sized", and "cluster-sized" 
ACDM halos. The radial dependence of the slope is significantly different 
from both the NFW and M99 profiles, and is better approximated by the p a 
profile described in eqs.(2) and (3). (b) Maximum inner asymptotic slope as 
a function of radius. Cusps as steep as p oc r -1 ' 5 are clearly ruled out by our 
simulations; (3q = 1 (as in the NFW profile) are still consistent with the data, 
although there is little evidence for convergence to a power law density profile 
at the center. 



Figure 2a shows the radial dependence of (5 in the simulated halos, and 
confirms that neither the NFW profile nor the M99 profile capture fully the 
structural diversity of ACDM halos. The profiles are clearly shallower near the 
center than the asymptotic value of (3o = 1.5 proposed by M99. At the same 
time, (5{r) deviates from the NFW prediction inside r„2 : density profiles become 
shallower less rapidly than expected from the NFW formula. Thus, NFW fits 
tend to underestimate systematically the density just inside r_2 (see Figure 
la). Fitting formulae with steeper cusps (such as the M99 profile) might indeed 
do better in that region; however, so will any other modification of the NFW 
formula where the logarithmic slope inside r_2 depends more weakly on radius 
than NFW's , regardless of (3q. 

Conclusive proof that the cusps cannot diverge as steeply as r -1,5 is provided 
by the total mass inside the innermost resolved radius, r conv (the minimum 
radius plotted for all profiles shown in this paper). This is because the mean 
density interior to any radius, p(r), together with the local density at that radius, 
p(r), provide a robust limit to the asymptotic central slope, /?o < /9 m ax(^) = 3(1— 
p{r)/p{r)) (under the plausible assumption that j3 is monotonic with radius). 
The radial dependence of /3 max is shown in figure 2b, and shows that, except 
for possibly one dwarf system, no simulated halo can have a cusp as steeply 
divergent as /?o = 1-5. 

Actually, there is no indication in the profiles for a well defined value for 
Po, as profiles keep getting shallower down to the smallest resolved radius. A 
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simple power law approximates the radial dependence of (3(r) better than the 
NFW or M99 profiles; 

P a (r) = -d\np/d\nr = 2(r/r- 2 ) a , (2) 
which corresponds to a density profile of the form, 

HPa/ P -2) = (-2/a)[(r/r^r-l}. (3) 

The thick dot-dashed curves in Figure 2 show that eq.2 (with a ~ 0.17) does 

indeed reproduce fairly well the radial dependence of (3{r) and (3 max {r) in sim- 
ulated halos. Furthermore, adjusting the parameter a allows the profile to be 
tailored to each individual halo, resulting in much improved fits over those ob- 
tained with the NFW or M99 formulae. The best-fit values of a (in the range 
0.1 - 0.2) show no obvious dependence on halo mass: the average a is 0.172 and 
the dispersion about the mean is 0.032 for the nineteen halos in our series (see 
Navarro et al 2003 for further details). 

3. CDM halos and LSB rotation curves 

A number of authors have reported disagreements between the shape of rotation 
curves of low surface brightness (LSB) galaxies and the circular velocity profiles 
implied by fitting formulae such as that proposed by NFW (Moore 1994, Flores 
k Primack 1994, McGaugh & de Blok 1998, de Blok et al 2001). Many of these 
galaxies are better fit by circular velocity curves arising from density profiles 
with a well defined constant density "core" rather than the cuspy ones inferred 
from simulations, a result that has prompted calls for a radical revision of the 
CDM paradigm on small scales (see, e.g., Spergel & Steinhardt 2000). However, 
before accepting the need for radical modifications to CDM it is important to 
note a couple of caveats that apply to the LSB rotation curve problem. 

• Strictly speaking, the observational disagreement is with the fitting for- 
mulae, rather than with the actual structure of simulated CDM halos. As 
noted in the previous section, there are small but systematic differences 
between them, so it is important to confirm that the disagreement persists 
when LSB rotation curves are contrasted directly with simulations. 

• It must be emphasized that the rotation curve problem arises when com- 
paring rotation speeds of LSB disks to spherically-averaged circular ve- 
locities of dark matter halos. Given that CDM halos are expected to be 
significantly non-spherical, some differences between the two are to be ex- 
pected. It is therefore important to use the full 3D structure of CDM 
halos to make predictions regarding the rotation curves of gaseous disks 
that may be compared directly to observation. 

Figure 3a illustrates the LSB rotation curve problem highlighted above. 
This figure shows the rotation curves of four LSB galaxies (points with error 
bars) selected from the sample of McGaugh et al (2001). The data points have 
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Figure 3. (a) Rotation curves of four LSB galaxies chosen to illustrate 
their various shapes, as measured by 7. The larger 7 the sharper the turn 
of the rotation curve from rising to flat, and the poorer the agreement with 
NFW profiles (dashed lines). The V c profile of simulated halos (solid lines) 
match reasonably well systems with 7 <J 2 but cannot account for those 
with 7 > 2.(b) Top panel shows the distribution of the "shape" parameter 
7 corresponding to all rotation curves in the samples of de Blok et al (2001, 
2002) and of Swaters et al (2003). Bottom panel compares the distribution 
of "observed" 7 with that derived from fits to halo circular velocities. This 
shows that it is possible to account for rotation curves with 7 < 2 but that 
curves with 7 J> 2.5 are inconsistent with halo V c profiles. 



been fitted using a simple formula, V vo t(r) = Vq(1 + (r/ro) -7 ) -1 / 7 (Courteau 
1997). Here Vq and ro are dimensional scaling parameters, whereas 7 is a di- 
mensionless parameter that characterizes the shape of the rotation curve. This 
three-parameter formula provides excellent fits to all galaxies, as witnessed by 
the quality of the (solid line) fits shown in figure 3a. 

To emphasize visually the shape discrepancy, the rotation curves in figure 
3a have been scaled to ro.3 and V0.3: the radius and velocity, respectively, where 
the slope of the curve is dlnVrot/dlnr = 0.3. The four galaxies have different 
values of 7, and have been chosen to illustrate to whole range of rotation curve 
shapes in the McGaugh et al and Swaters et al datasets. Roughly two-thirds of 
their LSBs have 0.5 < 7 < 2; the rest have 2 < 7 < 5 (see top panel of Fig.3b). 

The dashed lines in figure 3 show the V c profile of an NFW halo, which 
is fixed in these scaled units. Clearly, as 7 increases, the difference in shape 
between the NFW profile and the rotation curve becomes more acute. Galaxies 
with 7^2 are roughly consistent with NFW, whereas 7 > 2 rotation curves are 
clearly inconsistent. Can the deviations from NFW discussed in figures 1 and 2 
be responsible for this discrepancy? 

To address this question we have fitted all simulated halos with the same 
(Vo,ro,7) formula used for the observational data. Each halo is fitted at vari- 
ous times during its late evolution so as not to impose any restrictions on the 
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dynamical state of the halo; temporary departures from equilibrium associated 
with recent accretion events are thus taken into account 'naturally" in this pro- 
cedure. The 7 distribution of halos peaks at about 7 ~ 0.6 and is restricted to 
values that rarely exceed ~ 2.5 (see bottom panel of Fig. 3b). Halos whose fitted 
7 best match each of the galaxies in Figure 3a are shown with faint solid lines 
in that figure. 

Full details of this exercise are presented in Hayashi et al (2003), but the 
main conclusion is that it is possible to find halos whose V c curves match all 
galaxies with 7^2, with scaling parameters tq and Vq in reasonable agreement 
with observed values. CDM halos are thus actually consistent with the majority 
of LSB rotation curve data. Indeed, the rotation curve problem may be thought 
of as restricted to the few galaxies with 7 > 2; these rotation curves rise and 
turn too sharply to be consistent with the V c profiles of simulated halos (see the 
bottom right panel of figure 3a). 

Do these galaxies rule out the presence of a cusp in the dark matter density 
profile? As noted above, before concluding so we must take into account possible 
systematic differences between rotation speed and circular velocity in gaseous 
disks embedded within realistic, triaxial halos. This is a complex issue that 
involves a number of parameters, such as the degree of triaxiality, the role of the 
disk's self-gravity, size, and orientation, etc. 

For the sake of simplicity, we have decided to begin addressing this issue 
by evolving a massless, isothermal gaseous disk at the center of mildly triaxial 
(1:0.9:0.8) halos (Hayashi et al, in preparation). The triaxiality in the potential 
induces a radially-dependent inclination which may fluctuate, in the case we 
consider here, by 5 to 10 degrees from one radius to another. The effect on 
rotation curves of departures from circular symmetry and from coplanar orbits 
derived from long-slit spectra is complex, but the shape of the rotation curve, in 
particular, is affected. On some projections, rotation curves appear to rise and 
turn rapidly, and they would be (erroneously) taken to imply the existence of 
a constant-density core in simple models that impose spherical symmetry. Fits 
to the rotation curves using the (Vb,ro,7) formula often have 7 > 2, consistent 
with galaxies where NFW profiles provide a particularly poor fit to the rotation 
curve data. This suggests that deviations from spherical symmetry in the mass 
structure of CDM halos might reconcile rotation curve shapes that seem to favor 
the presence of constant density cores with cusps in the dark matter density 
profiles. 



4. Concluding Remarks 

Although it appears possible to reconcile dark matter cusps with LSBs by ap- 
pealing to asphericity in CDM halos, it would be premature to argue that the 
problem has been fully solved. After all, it may come as no surprise that one is 
able to reproduce LSB rotation curves, given the number of extra "free" param- 
eters afforded by relaxing the assumption of spherical symmetry. It is therefore 
important to build a more compelling case for this interpretation of rotation 
curve data, so as to render it falsifiable. We are in the process of identifying 
corroborating trends that may be used either to support or to rule out this inter- 
pretation, as the case may be. In particular, we would like to characterize better 
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the 7 > 2 systems: are such rotation curves in general asymmetric? Is rotation 
around the minor axis expected? Can one verify the triaxial-halo interpretation 
in full velocity maps? In this sense, identifying a "make or break" prediction 
will be as important as the success of aspherical halos in reproducing the rich 
variety of shapes of LSB rotation curves. Only once this is accomplished shall 
we be able to conclude that LSB rotation curves do not preclude the presence of 
dark matter density cusps, freeing the CDM paradigm of one of its most vexing 
challenges on small scales. 
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